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Abstract 

We investigate the effect of SiC stacking on the 4H-SiC/Si02 interface, both in the presence and 
absence of O defects, which appear during thermal oxidation, via first principles calculations. It 
is known that 4H-SiC(0001) has two different surface types, depending on which of the two lattice 
sites, h or k, is at the surface [K. Arima et ai, Appl. Phys. Lett. 90 , 202106 (2007)]. We find 
interlayer states along the conduction band edge of SiC, whose location changes depending on the 
interface type, and thus too the effect of defects. When h sites are directly at the interface, O 
defects remove interfacial conduction band edge states. On the other hand, when k sites are at the 
interface, the conduction band edge is insensitive to the presence of O defects. These differences 
will impact on the operation of SiC devices because the most commonly used SiC based metal- 
oxide-semiconductor field-effect transistors rely on the electronic structure of the conduction band. 
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Silicon carbide (SiC) is a wide band-gap semiconductor material that can be used for 
high temperature and power electronic devices, in conditions where the traditionally used Si 
fails. These improved properties are due to the presence of strong Si-C bonds. Like Si, its 
native oxide is Si 02 , which can be grown via thermal oxidation, making it useful for metal- 
oxide-semiconductor held-effect transistors (MOSFETs), primarily of the n-channel type. 
However, its use in practical devices has been hampered by the low channel mobility of the 
SiC/Si02 interface compared with bulk SiGi*^. This has been partly attributed to a variety 
of C and O defects found on both sides of the interfacial region^^-. Traditionally research has 
focused on defects which contribute interface gap states, but advances in interface growth 
techniques mean that these gap states no longer make a signihcant contribution to reductions 
in mobility^'-. Further advances require new information. 

SiC has many different polytypes, the most commonly used of which is 4H-SiC. This 
consists of four repeated SiC bilayers, with two inequivalent lattice sites for the Si and C 
atoms. Based on their local symmetries, these are often labelled h for hexagonal and k for 
quasi-cubic, as shown in Fig. [TJ At h sites the connected SiC bilayers have Si-C bonds in 
opposite directions, whilst for k sites these Si-C bonds run parallel. The local structure 
around the SiC/Si02 interface will change depending on which of these sites is at the SiC 
surfacel^iii. For convenience these will be referred to as h and k type interfaces based on 
the site which the surface Si atom occupies. Previous work looking at bulk SiC reported 
interlayer floating states at the conduction band edge (CBE)^J^, but have not reported the 
effect of defects at the SiC/Si02 interface. Floatings states appear between subsequent SiC 
bilayers with the same orientation. 

Here we investigate how the energies and local density of states (LDOS) of the inter¬ 
face changes with interface type, and when O atoms are subsequently introduced. This is 
important for real devices, because currently the surface SiC bilayer cannot be controlled 
during the thermal oxidation process. Due to differences in local structure, floating states 
are observed at the interface for h type, but from the second bilayer down for k type. More 
interestingly, O defects behave differently in the two interfaces, removing the floating states 
at the h type whilst leaving the k type unchanged. The variation of the floating states at 
the CBE is expected to affect the carrier mobility at the interface because the operation of 
n-channel MOSFETs relies on the electronic structure at the CBE. 

Calculations are carried out using RSPACEl^, which performs density functional the- 


2 


oryisj^ calculations using a real-space finite-difference methodic, employing a time-saving 
double-grid technique^^^— . The electron-ion interactions are described using the projector- 
augmented wave method^i for the C, Si and O atoms, and a norm-conserving pseudopo- 
tential^^i^ for H. The exchange-correlation interaction is treated within the local density 
approximation^^. 

The interface is modelled using 6 bilayers of 4H-SiC and 9 A thick of /9-tridymite Si02, 
with the interface made at the SiC(OOOl) surface, saturating all the Si dangling bonds with O. 
Both the top and bottom layers of the model are H terminated. The interface is constructed 
using experimental values for the Si-C and Si-0 bond lengths. The lateral size of the cell is 
SiC(OOOl) 3 X 3\/3 with the z direction perpendicular to the SiC surface. The cell is periodic 
in all directions with an 11 A vacuum gap in the z direction. Calculations are performed 
with a coarse grid spacing of 0.16 A and using the T and M k-points of SiC(OOOl) 3 x 3\/3. 
The bottom H atoms and SiC bilayer are held hxed and structural optimization is performed 
on all remaining atoms to reduce the forces between atoms. A variety of known O defects^'^- 
are introduced into the SiC. Full details on interface atomic structures including defects will 
be described later. 

Formation energies for defects are calculated according to 

^form T Po ^nO (f) 

where Eno is the energy of the interface containing n excess O atoms and po is the chemical 
potential of oxygen, calculated from the energy of an O 2 molecule. The energies of single 
O defects are compared to the clean interface, whilst double O defects are compared to the 
lowest energy single O defect. 

Formation energies can reveal which defect structures are more likely to form, giving 
an idea of their relative abundance at the interface, but do not give an insight how their 
properties differ from the clean interface. Thus we also examine the LDOS for both types 
of interface, with and without defects. The LDOS is calculated according to 

p{z,E) =Y, [ \^iA^,y,z)fdxdy X (2) 

i,k 

where Si^k are the eigenvalues of the wavefunction, with indexes of i and k for the eigenstate 
and the k-point respectively. N = is the normalisation factor, with a as the smearing 
factor, here set to 13.5 eV“^. 
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Before looking at defects, it is useful to compare the clean interface. We find that the 
h type interface is lower in energy than the k type by 0.27 eV per 3 x 3\/3 region. The 
LDOS for these two structures, as shown in Fig. O is mostly the same, except for along 
the CBE, where several oval-like features can be seen, whose location changes with interface 
type. For h type, one of these states is at the interface, whereas for k type the hrst state 
does not appear until the second bilayer. This changes the band gap in the vicinity of the 
interface. These CBE states are the previously noted floating states^^ii^, with the difference 
in location a reflection of the local structure at the interface. To aid with future discussions 
on this matter, we introduce channel depth, a measure of the number of subsequent SiC 
bilayers with the same orientation. Channel depth is essentially the same concept as the 
channel length discussed in Ref. Il3|, but the procedure for counting the length of a channel is 
different. We use a slightly different term because our counting procedure is more intuitive 
for describing the interface. For bulk 4H-SiC the longest channel depth is two. At the h 
type interface the channel depth is also two, but for the k type interface the depth is only 
one. A minimum channel depth of two is required to observe floating states. This distinction 
between interface types is important, especially given the position of interfacial defects. 

We then introduce O atoms into the interface as shown in Fig. |3l In our previous study^, 
we investigated the total energies of all possible O atom sites and obtained the most stable 
conhgurations. Here, the lowest energy structures which appear during one oxidation cycle 
are examined. For a single excess O this is Oij^ [Fig. [3][a)], an interstitial site at the 
interface. For two excess O atoms this is Osub+Oif^ [Fig. |3][b)], which combines subsurface 
and interstitial sites. Post CO emission this is VcO^ [Fig. [3][c)], which involves O atoms 
passivating the dangling bonds created by the C vacancy. 

The formation energies of defects are calculated for both interface types, with energy 
orderings agreeing with prior calculations^. Formation energies for Oj/ are 0.17 and 0.25 eV 
for the h and k type interfaces respectively. For Ogub+Oif these formation energies are 0.60 
and 0.83 eV. In most cases formation energies are greater in the k type interface, despite the 
clean interface being higher in energy. This means k type interfaces would have a slightly 
higher defect density, but as the following LDOS show, this is of little consequence because 
its electronic structure is largely insensitive to O defects. 

In addition to their energies, we also examine the LDOS for each defect structure. Firstly 
Oj/, as shown in Fig. 01 where the most signihcant changes occur at the CBE, and only for 
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the h type interface. The hrst floating state is removed, shifting the position of the CBE 
and increasing the band gap. The presence of an extra O atom has a channel blocking 
effect, reducing the channel depth at the interface, and removing the floating state. The k 
type interface shows no signihcant changes at the CBE because there is no floating state 
to remove. It is also useful to study a single subsurface site, Oauh [Fig- El[d)], to see how O 
location influences the channel blocking effect. The same changes at the CBE are observed 
for Osuh (not presented here), meaning that defects located within a bilayer or between two 
bilayers both have a channel blocking effect. On the other hand, changes at the valence 
band edge (VBE) between the clean interface and Ojj are the same for both interface types. 
A defect state is observed at the VBE, but the state is not distinct from the bulk SiC states. 
This defect state can be attributed to a C-0 tt bond, as observed in the partial DOS for 
the highest occupied band (not presented here). Because the dissociation energy of a C-0 
bond at the SiC/Si02 interface is lower than that of a Si-C bond,"^ the VBE at the interface 
shifts up in energy slightly. 

Secondly, Ogub+Oif as shown in Fig. [5l Changes at the CBE are in line with earlier 
results, with the interface floating state removed for h type, but k type mostly unaffected. 
This supports the role of O defects in channel blocking. Changes at the VBE are more 
noticeable than previous results, with a defect state which is distinct from the bulk SiC. 
Examining the partial DOS for the highest occupied bands shows that this defect state 
again arises from C-0 vr bonds. The important difference between Osub+Oif and 0^/ is that 
it contains an O-C-0 bond, which the latter does not. As such, the states associated with 
Osnfe+Oj/ are even higher in energy, making them distinct from the VBE. 

Finally, the LDOS for the Vc02 structure, as shown in Fig. [6l which looks similar to 
the clean interface. Once again, at the CBE the floating state at the h type interface is 
removed due to channel blocking by the defect, with no change for the k type. In addition, 
no changes are observed at the VBE for both interface types. There are no defect states 
because there are no C-0 bonds present. 
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We have also examined some of the C-related defect structures reported in Refs, 
and found defect gap states which appear close to the CBE, partially obscuring changes 
to the floating states. However, if results for C-related defects are included our overall 
conclusions do not change. 

These results are summarized in Fig. [TJ with the primary focus on the h type interface. 
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since it undergoes greater changes in the presence of defects. The clean h type interface is 
represented by (a). When defects are introduced, the band diagram is modihed in two main 
ways, as shown in (b), (c), and (d) for the O*/, Ogub+^if-, and Vc 02 defects respectively. 
Firstly, states at the CBE are removed, opening the band gap near the interface. This occurs 
due to channel blocking by O defects, and occurs in all cases reported here. The gap at 
the CBE extends over the hrst two bilayers because the state that would appear there is 
removed. Secondly a defect state can appear near the VBE, reducing the band gap. For 
single O atoms, such as O*/, the band edge at the interface shifts up in energy slightly, 
as shown in (b). The defect state is only distinct from the bulk SiC states in the case of 
Osn6+Oj/, as shown in (c). 

A similar gap is observed at the CBE of the k type interface, whether defects are present 
or not, but with a smaller extent. The gap only extends over the hrst bilayer, a rehection 
of the channel depth at the interface. If the extent of this gap is ignored, both the clean 
interface and the Vc02 defect can be depicted by Fig. [7|^d), with the Oj/ and Ogub+Oif 
defects represented by (b) and (c) respectively. However, diagram (a) never appears for the 
k type interface. 

In summary, we hnd that interface type determines whether CBE hoating states at the 
interface are affected by the presence of O defects, as with h type, or not, as with k type. 
Floating states are observed between adjacent SiC bilayers with parallel Si-C bonds and 
are found from the hrst bilayer down for the h type interface and from the second bilayer 
down for the k type. When O defects are introduced at the interface they have a channel 
blocking ehect, removing the hoating state at the h type interface, with little ehect on the 
k type interface, where there are no states to remove. Defect states are also found at the 
VBE, arising from the C-0 tt bonds formed by the O defects, but are only distinct from 
the bulk SiC states in the case of the Ogub+Oif defect. These states do not vary between 
the h and k type interfaces. Overall the CBE of the h type interface is more ahected by 
defects, so it would be expected that hoating states will inhuence the electronic properties 
of n-channel MOSFET, such as carrier mobility. In the future we want to do more complete 
calculations for all of these structures, to get a fuller picture of how the properties of the 
interface change with defect type. 
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FIG. 1. (Color online) SiC/Si02 interface with (a) h and (b) k lattice sites at the interface. Black 
lines indicate the bounding box of the calculation cell. Highlighted Si-C pairs show examples of h 
and k sites. Si atoms are beige, C grey, O red and H white. 
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FIG. 2. (Color online) LDOS for the h (left) and k (right) type interface with no defects. The 
horizontal axis is energy relative to the Fermi energy, Ei?, defined as the middle of the band gap. 
The vertical axis is the height of the model. For clarity, structural models are provided to the 
right of each LDOS. Contours correspond to density, from blue at the lowest to red at the highest, 
doubling at each contour line, from a minimum of 5.3x10“^ states eV“^A“^. 
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FIG. 3. (Color online) Structural models for the various O defects in the h type interface, (a) Oj/ 
from above, (b) Osub+Oif from both above and the side, (c) Vc 02 from above, (d) Osub from the 
side. Views from above look down on the SiC(OOOl) surface, with the Si02 above the connecting 
region removed for clarity. Black boxes indicate the bounding box of the calculation cell. Side 
views looking down the [1100] direction are to clarify the position of subsurface atoms. Arrows 
indicate defect positions. 
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FIG. 4. (Color online) LDOS for the Ojj defect in the h (left) and k (right) type interface. Notation 
is the same as Fig. [2l These LDOS focus on the region around the SiC band gap. 
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FIG. 5. (Color online) LDOS for the Ogub+Oif defect in h (left) and k (right) type interface. 
Notation is the same as Fig. [2j These LDOS focus on the region around the SiC band gap. 
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FIG. 6. (Color online) LDOS for the Vc02 defect in h (left) and k (right) type interface. Notation 
is the same as Fig. [2l These LDOS focus on the region around the SiC band gap. 
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FIG. 7. Schematic band diagrams for the h type SiC/Si02 interface, highlighting important changes 
caused by defects, (a) Clean interface, (b) Single O defect, (c) Double O defect (d) Post CO 
emission. 
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